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ABSTRACT 

Aims. We performed an integral field spectroscopic study for the HII galaxy IIZw70 to investigate the interplay between its ionized 
interstellar medium (ISM) and the massive star formation. 

Methods. Observations were taken in the optical spectral range from /13700 A - 6800 A with the Potsdam Multi-Aperture 
Spectrophotometer (PMAS) attached to the 3.5 m telescope at Calar Alto Observatory. We created and analyzed maps of spatially dis- 
tributed emission-lines (at different stages of excitation), continuum emission, and properties of the ionized ISM (ionization structure 
indicators, physical-chemical conditions, dust extinction, kinematics). We investigated the relation of these properties to the spatial 
distribution and evolutionary stage of the massive stars. 

Results. For the first time we have detected Wolf-Rayet (WR) stars in this galaxy. The peak of the ionized gas emission coincides 
with both the location of the maximum of the stellar continuum emission and the WR bump. The region of the galaxy with lower 
dust extinction corresponds to the region that shows the lowest values of velocity dispersion and radial velocity. The overall picture 
suggests that the ISM of this region is being disrupted via photoionization and stellar winds, leading to a spatial decoupling between 
gas+stars and dust clouds. The bulk of dust appears to be located at the boundaries of the region occupied by the probable ionizing 
cluster. We also found that this region is associated to the nebular emission in HeIL14686 and to the intensity maximum of most 
emission lines. This indicates that the hard ionizing radiation responsible for the HeIL14686 nebular emission can be related to the 
youngest stars. Within ~ 0.4 x 0.3 kpc 2 in the central burst, we derived oxygen abundances using direct determinations of r e [OIII]. 
We found abundances in the range 12+log(0/H)= 7.65 - 8.05, yielding an error-weighted mean of 12+log(0/H)=7.86 ± 0.05 that has 
been taken as the representative oxygen abundance for IIZw70. 

Key words. ISM: HII regions - ISM: ionization structure - Massive stars - Galaxies: abundances - Galaxies: dwarf - Galaxies: 
individual (IIZw70) 



1. Introduction 

The HII galaxies are gas-rich, dwarf systems that have expe- 
rienced intense recent or ongoing violent star formation (SF). 
These objects were identified for the first time by Haro (1956) 
and Zwicky (1964), and are characterized by the presence of gi- 
ant HII regions that dominate their observable properties at op- 
tical wavelengths. In particular, their integrated spectra are very 
similar to those of giant extragalactic HII regions in spiral galax- 
ies (Sargent & Searle 1970). The analysis of their emission line 
spectrum indicates that HII galaxies are metal poor objects (1/40 
Z < Z < 1/3 Z Q ; e.g. Kunth & Sargent 1983). 

HII galaxies are ideal laboratories for probing the interplay 
between massive SF and the ISM in low metallicity environ- 
ments. The massive SF process in dwarf galaxies has effects on 
the properties of the surrounding ISM. A large number of mas- 
sive stars (between 10 4 - 10 6 solar masses of gas are transformed 
into stars in HII galaxies) are formed almost simultaneously 
within relatively small volumes (e.g. Melnick 1992). Millions of 
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years after the onset of the burst, the most massive stars begin to 
explode as supernovae, creating violent, short-lived injections of 
kinetic energy and metallic elements into the ISM. The disrup- 
tion of the ISM may significantly affect the spatial distribution 
of gas and dust particles in the regions close to the massive star 
cluster, and also determine the way in which new metals ejected 
by massive stars are mixed with the original gas from which the 
stars formed (Kunth & Ostlin 2000). 

High spatial resolution imaging has revealed that in many 
HII galaxies the ionized material presents a complex structure: 
star clusters in the main body (with a non-uniform distribution 
of SF knots, ensembles of star clusters, or individual super stel- 
lar clusters) of the galaxy and external gas components out- 
side the young stellar clusters (e.g. Martin & Kennicutt 1995, 
Vflchez & Iglesias-Paramo 1998, Papaderos et al. 2002). Thus, 
a two-dimensional analysis of the ionized material in HII galax- 
ies yields spatially resolved information on properties of the ion- 
ized gas. The study of the distribution of these properties is an 
important issue for our understanding of the interplay between 
the massive stellar population and the ISM. 
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Fig. 1. Image of IIZW70 in the r (left) and u (right) bands from 
Sloan Digital Sky Survey. The blue box mosaic represents the 
field used in this work (left), which is located approximately in 
the central parts of the galaxy; the green box marks the likely 
ionizing clusters location, as indicated with a box on the u image 
(right). North is up, east is left. The spatial scale is shown in each 
panel, where 1 kpc ~ 11" assuming a distance to the galaxy of 
18.1 Mpc (see table 1). 



However, up to now only a few two-dimensional spectro- 
scopic studies of HII galaxies have been performed. Recent work 
by Izotov et al. (2006) presents two-dimensional spectroscopy of 
the extremely metal-deficient HII galaxy SBS 0335-052E. They 
found a small gradient of the electron temperature T e and oxy- 
gen abundance, and the presence of an ionized gas outflow in the 
perpendicular direction to the galaxy disk. Cairos et al. (2002), 
using two-dimensional spectroscopy, obtained the ionized gas 
velocity field in the central part of Mrk 370. 

In this paper, we present a two-dimensional spectroscopic 
study for the HII galaxy IIZw70. This galaxy is number 70 
in Zwicky's second list of "Compact Galaxies and Compact 
Parts of Galaxies, Eruptive and Post-eruptive Galaxies"(Zwicky 
& Zwicky 1971). IIZw70 has been classified as an HII galaxy 
(Sargent 1970). The basic data of IIZw70 are shown in table [1] 
This object is interacting with its nearby companion IIZw71 at a 
projected distance of 23 kpc (assuming a distance to the IIZw70- 
71 system of 18.1 Mpc; see table 1). The ongoing interaction is 
indicated by interferometric HI studies (Balkowski et al. 1978, 
Cox et al. 2001) which revealed a gaseous streamer connecting 
IIZw70 with the polar ring galaxy candidate IIZw7 1 (Whitmore 
et al. 1990, Reshetnikov & Combes 1994). Cairos et al. (2001) 
show that IIZw70 presents very elongated outer isophotes and a 
very blue (U-B = -0.89) nuclear starburst in what seems to be 
an edge-on disk. They also find that the starburst activity is con- 
centrated in the optical center of the galaxy and is surrounded by 
faint gaseous emission with quite a distorted morphology. 

In figure Q] we show the images of IIZw70 in the r (left im- 
age) and u (right image) bands extracted from SDSSQ. Over the r 
filter image, we marked the field of view used in this work (blue 
box mosaic). 

In this work, we investigate the relation between massive SF 
process and the ionized ISM in IIZw70 using integral field spec- 
troscopy (IFS). In order to do this we created and analyzed maps 
of spatially distributed emission-lines, continuum emission and 
physical-chemical properties of the ionized gas (electron tem- 
perature and density, gaseous metal abundances, dust extinction, 
excitation). We present a discussion on the significance of the 
chemical abundance variation found here. We also present a brief 
study of the kinematic properties of the ionized gas. The rela- 
tionship between the physical-chemical and kinematic proper- 
ties spatial distribution is discussed. 



1 The Sloan Digital Sky Survey website is http://www.sdss.org 



Parameter 


Value 


Name 


IIZw70 


Other designations 


UGC 9560, Mrk 829 


R.A. (J2000.0) 


14h 50m 56.5s 


DEC. (J2000.0) 


+35d34' 18" 


redshift 


0.004 


M^mag) 


-16.3 


m^mag) 


15.1 


D c (Mpc) 


18.1 




0.34 x 10" 




1/5 



" Absolute magnitude in B from Cairos (2001) 

b Apparent magnitude in B from Deeg et al. (1997) 

c Distance to the galaxy from Cox et al. (2001) 

d Neutral hydrogen mass from Thuan & Martin (1981) 

e Metallicity from Kobulnicky & Skillman (1996) 



Table 2. Log of the observations. 



Pointing 


Exptime 


Grating 


Spec, range 


Date 


IIZw70 


(s) 




(A) 




IIZw70-NE 


3x750 


V300 


3900-7049 


05/08/2005 


IIZw70-SW 


3x750 


V300 


3900-7049 


05/08/2005 


IIZw70-NE 


3x1200 


VI 200 


4330-5020 


06/08/2005 


IIZw70-SW 


3x1200 


VI 200 


4330-5020 


06/08/2005 


IIZw70-NE 


3x1200 


VI 200 


3610-4379 


04/06/2006 


IIZw70-SW 


3x1200 


VI 200 


3610-4379 


04/06/2006 



In the following section we describe the observations and 
data reduction. In Sect. 3 the results are presented. We discuss 
our results in Sect.4. 

2. Data 

2.1. Observations 

IFS allows us to simultaneously collect spectra of many different 
regions of an extended object, under identical instrumental and 
atmospheric conditions. This is not only a more efficient way of 
observing, but it also guarantees the homogeneity of the dataset. 

The spectra of IIZw70 were taken with PMAS that is an 
integral field unit (IFU) developed at the Astrophysikalisches 
Institut Potsdam (Roth et al. 2005). It is currently installed at 
the Calar Alto Observatory 3.5 m Telescope. Two pointings of 
the galaxy (IIZw70-NE and IIZw70-SW) were taken in order to 
cover the central regions and different regions of interest of the 
galaxy (e.g. the tails). 

The log of observations is given in table [2] The V300 and 
VI 200 gratings present a dispersion of 1.6 A/pixel and 0.36 
A /pixel, respectively. The PMAS spectrograph is equipped with 
256 fibers coupled to a 16 x 16 lens array. Each fiber has a 
spatial sampling of 1" x 1" on the sky resulting in a field of 
view of 16" x 16". Calibration images were obtained follow- 
ing the science exposures and consisted of emission line lamps 
spectra (HgNe), and spectra of a continuum lamp needed to lo- 
cate the 256 individual spectra on the CCD. Observations of the 
spectrophotometric standard stars BD +28°421 1 and Hz44 were 
obtained during the observing nights for flux calibration. The 
nights were photometric and the typical seeing during the obser- 
vations was 1". 
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2.2. Data reduction 



We reduced the data using the software R3D (Sanchez 2006). 
Different exposures taken at the same pointing were combined 
using IRAF0 tasks. The expected locations of the spectra were 
traced on a continuum-lamp exposure obtained before each tar- 
get exposure. After bias subtraction, we extracted the target 
spectra by adding the signal from the 5 pixels around the central 
traced pixel (that is the total object spectrum width). With expo- 
sures of Hg and Ne lamps obtained immediately after the science 
exposures, the spectra were wavelength calibrated. We checked 
the accuracy of the wavelength calibration using sky emission 
lines, and find standard deviations of 0.5 A and 0.3 A for 
the V300 and VI 200 gratings, respectively. The continuum-lamp 
exposure was also used to determine the response of the instru- 
ment for each fiber and wavelength (the so-called fiat-spectra). 
We used these flat-spectra in order to homogenize the response 
of all the fibers. The average nearby background spectrum, ob- 
tained moving the IFU off-target, was subtracted from the tar- 
get spectra. For the standard star observations we co-added the 
spectra of the central fibers and compared the one-dimensional 
standard star spectrum with table values to create a sensitivity 
function. The spectra were flux calibrated using IRAF. 

The reduced spectra were contained in two data cubes, each 
one corresponding to a different PMAS pointing. The two data 
cubes were mosaiced into a single frame (see the blue box mo- 
saic in figure Q]). We corrected the mosaic data cubes for the ef- 
fect of differential atmospheric refraction using the R3D pack- 
age (Sanchez 2006). We recentred and scaled the mosaic data 
cube taken in 2006 with respect to the mosaic data cubes ob- 
tained in 2005 making use of the Hy emission line, common 
to the wavelength ranges of both gratings. Further manipulation 
of the cubes, such as the production of emission-line maps was 
performed using E3D (Sanchez 2004). 

2.3. Line intensities 



og f(Ha) 




10 15 20 

A RA (arcsec) 



Ha/H(S 




5 10 15 20 

A RA (orcsec) 



Fig. 2. Top image: map of the log flux(Ha) in each pixel in units 
of 10~ 15 erg s" 1 cm -2 ; Ha fluxes with relative errors > 10% are 
not shown. We overplot the exact location of the green box from 
figure 1 for orientation. Bottom image: Ha/H/3 ratio spatial dis- 
tribution. The lowest values tend to be concentrated near to the 
brightest region of the galaxy. Only HafH/3 ratios with relative 
errors < 30% are shown. North is up, east is left. Relative right 
ascension (RA) and declination (DEC) are shown in arcsec. 



Using the software FIT3D (Sanchez et al. 2007), we fitted line 
profiles on the extracted one-dimensional spectra in order to de- 
rive the integrated flux of each emission line. A single gaussian 
was fitted to each emission line, using a polynomial function to 
characterize the continuum. The software FIT3D allows defini- 
tion of emission line systems, i.e., a kinematically coupled set 
of emission lines with the same width. This was essential for 
accurate deblending of the lines, when necessary. 

For each fiber spectrum we derived its corresponding redden- 
ing coefficient, C(H/3), using the value of the Balmer decrement 
derived from Ha/H/3, as compared to the theoretical value ex- 
pected for case B recombination from Storey & Hummer ( 1 995), 
and applying the extinction law given by Whitford (1958). Thus, 
the fluxes of the emission lines for each fiber were corrected for 
extinction using their corresponding C(H/3) value. 

Equivalent widths (EWs) of Balmer emission lines are so 
high [e.g. EW(H/3) ~ 40-120 A] that the effect of the underly- 
ing stellar population in these lines appears not to be important 
[underlying absorption-lines EWs present typical values of 1-2 
A ; see e.g. McCall et al. (1985)]. In addition, we carried out an 
eye-inspection, fiber-by-fiber, of the H/3 and Hy emission lines, 
and we did not find any stellar absorption features (e.g. wings of 
absorption lines) underlying these emission lines. 



2 IRAF is distributed by the National Optical Astronomy 
Observatories. 



We calculated the error in the line fluxes, 07, from the ex- 
pression cr, = o>N 1/2 [l+EW/NA] 1/2 (Catellanos 2000) where a c 
represents a standard deviation in a box centred close to the mea- 
sured emission line, N is the number of pixels used in the mea- 
surement of the line flux, EW is the equivalent width of the line, 
and A is the wavelength dispersion in A/pixel. This expression 
takes into account the error in the continuum and the photon 
counts statistics of the emission line. The error measurements 
were performed on the extracted one-dimensional spectra. 

Ho- emission line map (continuum subtracted) and the spatial 
distribution of Ha/H/3 are shown in figure [2] The dust distribu- 
tion seems to be non uniform across the galaxy. We would like to 
stress that in the following sections, all emission line fluxes, used 
to derive the ionization structure and physical-chemical param- 
eters, are corrected for extinction using the corresponding value 
ofC(Hy6)foreach fiber. 

3. Results 

3.1. Galaxy structure 

In figure [3] we show the continuum emission near H/3 and Ha 
maps, together with the corresponding H/3 and Ha emission 
line intensity contours overplotted. The galaxy's morphology ap- 
pears very similar in emission lines and in continuum emission, 
i.e. the continuum and emission line maps (as seen from the con- 
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Fig. 3. Maps in the continuum near Hj3 (left) and Ha (right) 
emission lines; isocontours of the corresponding emission line 
fluxes are shown overplotted. Maps of the continua are repre- 
sentative of the stellar emission, free of contamination from the 
gaseous emission lines. Isocontours and maps are displayed in 
logarithmic scale, with north up and east to the left. The image 
size is ~ 24" x 23". The pixel size is ~ 0.6". 




Fig. 4. The zone occupied by the probable brightest clusters 
(blue rectangular box with physical dimension as green box in 
figure[T]( right)), detections of the WR bump (green crosses) and 
nebular Hell/14686 emission (blue crosses) on Ha flux map (see 
text for details). North is up, east is left. The size of the image is 
25" x 21" and the pixel size is 1". The spatial scale is shown, 
where 1 kpc ~ 11" assuming a distance to the galaxy of 18.1 
Mpc (see table 1). 

tours in figure |3]l show a central region with ionized gas and a 
faint extension towards the south-west. However we note that in 
the emission line contours we can resolve clearly two filamen- 
tary structures, with a tail-like shape, at the faintest extensions 
of the galaxy which are not present in the continuum maps. 

The r filter image in figure Q] shows a red extension already 
found by Cairos (2001) and Kong et al. (2003). Such a red ex- 
tension could indicate the likely existence of an intermediate- 
old underlying stellar population. Kong et al. (2003) studied the 
stellar content of a sample of 73 HII galaxies by analysing their 
continuum spectra and absorption features. They showed that in 
most of HII galaxies (the galaxy IIZw70 is among them), stars 
older than 1 Gyr exist. The average colours of the low-surface- 
brightness component (B-V = 0.47, V-I = 0.89) are also consis- 
tent with an evolved stellar background with an age of 1-3 Gyr 
forIIZw70, as shown by Cairos (2001). 

The u band (right image in figure [TJ is expected to trace the 
youngest stellar populations. The box in this image coincides 
with the u band maximum and encloses the region occupied by 
the probable brightest young clusters. We also show the exact 
location of this box in the r image of IIZw70 as a reference (left 
image). We verified that the location of the probable youngest 
stellar population, as seen in the u image, coincides with the 



central bright SF knot and with the zone of maximum in the 
continuum emission (see figures[3]and 0). 

3.2. Ionization structure and excitation sources 



I: g [SII]/Ha I r [Ml] Ha 




Fig. 5. Emission line ratio maps in logarithmic scale: 
[Sn]^6717.31/Ha (top left panel), [NII]^6584/Ha (top right 
panel) and [OIII]/l5007/Hy8 (bottom panel). Only line-intensity 
ratios with relative errors < 30% are represented. Isocontours of 
the Ha emission line flux are shown overplotted as reference. 
North is up, east is left. Relative right ascension (RA) and decli- 
nation (DEC) are shown in arcsec. 

The line ratio maps for the usual diagnostic lines (e.g. 
Baldwin et al. 1981; Osterbrock & Ferland 2006) can be used 
to investigate the spatial distribution of the ionization struc- 
ture of the ionized nebulae. Figure [5] shows the line ratio maps 
(in logarithmic scale) for [NII]^6584/Hff,[SII]^6717,31/Ha, 
[OIII]/15007/H/3. Ha emission line flux contours are plotted over 
all maps. High excitation values correspond to low values of the 
first two ratios and to high values of the last one. 




-1.5 -1.0 -0.5 0.0 -1.5 -1.0 -0.5 0.0 

log [Sll]6717,31/Hn log [Nll]6584/Ha 

Fig. 6. Left panel shows the relation between [OIII]/15007/Hj8 
and [SII]/l6717.31/Hff, and right panel presents the relation be- 
tween [OIIIM5007/H/? and [NII]/16584/Ha in each fiber. The 
solid lines are adapted from Osterbrock & Ferland (2006) (see 
the text for details). 
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It can be seen that [NII]/}6584/Ha as well as 
[SII]/l6717,31/Ha increases in zones where the Ha inten- 
sity decreases; [OIII]/15007/H/? presents the opposite trend. This 
behaviour towards lower excitation when the intensity is smaller 
can be related to the increasing distance from the young stellar 
clusters, with the corresponding decrease in the local values of 
the ionization parameter, as it is typical in extended HII regions 
(see McCall et al. 1985) 

Figure [6] shows the classical BPT diagnostic diagrams 
(Baldwin et al. 198 1) using our own data. Only fibers that present 
diagnostic line-intensity ratios with relative errors < 30% are 
plotted. The solid lines (adapted from Osterbrock & Ferland 
2006) show the locus of separation between regions dominated 
by photoionization (Hll-regions like; left of the line) and regions 
dominated by shocks. Our data show that the line ratios for most 
positions in the galaxy are located in the general locus of HII- 
region like objects. 

There are a few fibers, located at the tail-like regions 
(see the [SII]i6717,31/Ha map in figure 0, with values of 
[SII]/t6717,31/Ha somewhat higher that are found outside the 
HH-region zone. These values are consistent with predictions 
from shock models (Dopita & Sutherland 1995) suggesting the 
presence of shocks in the ionized gas outflow scenario proposed 
in section 3.5. Nonetheless, we must keep in mind that the line 
ratios of both axis plotted in figure[6]are sensitive to the metallic - 
ity, as well as, to the ionization parameter. Besides, as these line 
ratios can vary within an HII galaxy (i.e. across giant HII re- 
gions), thus BPT diagnostic diagrams on spatially resolved data 
of HII galaxies should be studied in depth though that is not the 
purpose here. 

3.3. Wolf-Rayet stellar population 

WR stars are the evolved descendants of massive O progenitor 
stars (Maeder & Conti 1994). Stellar evolution theory predicts 
few WR stars to form in a low-metallicity environment such 
as in IIZw70 (e.g. Meynet 1995). Therefore WR star detections 
in low metallicity HII galaxies have important implications for 
the stellar evolution theory. A common property of HII galaxies 
with WR stars (WR galaxies), provided that our understanding 
of stellar evolution is correct, is the existence of ongoing active 
or recent SF which has produced stars massive enough to evolve 
to the WR stage. This typically indicates ages of a few Myr for 
stars with initial masses M ini > 20M o (Maeder & Conti 1994). 

The presence of WR stars can be recognized via the WR 
bumps around /14650 A (blue bump) and /15808 A (red bump), 
which are generally a blend of Hell and several metal lines. The 
blending of several stellar and nebular emission lines around 
,14686 A and the complex spatial morphology in ground-based 
data make it challenging to disentangle the emission of stars and 
gas and to derive the WR content and the corresponding nebular 
properties. 

For the first time we detected WR stars in IIZw70. Summing 
the emission from 4 fibers (green crosses in figure|4]i we can ob- 
serve the blue and red bumps, - the latter marginally detected 
-, in the corresponding spectrum (see figure |7). This detection 
would indicate that WR bumps are being seen in an extended 
area (~ 4 arcsec 2 ). The use of the IFS minimises the difficulty 
in correlating spectra and image location, shown in other works 
(e.g. de Mello et al. 1998). The WR bumps intensity maximum 
coincides with the likely youngest clusters location (as shown 
in the u SDSS image; see figures Q] and [U and, within 1", also 
spatially coincides with the peak of the Ha emission. In addition 




I i I i I i LJI 

4000 4500 5000 5500 

Wavelength (A) 



Fig. 7. Spectrum showing the WR bumps obtained adding the 
emission from the 4 fibers (2" x 2") marked with a green cross 
in figure |4] Dashed line is the continuum fit. 



we verified that the maximum of the continuum emission adja- 
cent to H/3 (~ /14800 A) corresponds spatially to the WR bumps 
intensity maximum. From the spectrum shown in figure [7] we 
measured a value of EW(WRbump) of -2.2 + 0.6 A for the WR 
blue bump and a corresponding 100 x I(WRbump)/I(HjS) of 4.6 
+ 1.1 . 

An area of Hell/14686 narrow emission was found in 3 fibers 
located in the southwest corner of the box encompassing the 
probable central brightest clusters (see figure |4); the line is nar- 
row showing approximately the same width as the [OIII] lines, 
confirming its correct identification as a nebular line. The pres- 
ence of the nebular Hell/14686 line should indicate the existence 
of a very hard ionizing spectrum. 

3.4. Physical properties and abundance analysis 

We calculated the physical properties and ionic abundances of 
the ionized gas for IIZw70 following the 5-level atom FIVEL 
program (Shaw & Dufour 1994) available in the task IONIC of 
the STSDAS package. Total abundances were derived as we de- 
tail in the following. We calculated the final quoted errors in the 
derived quantities by error propagation taking into account er- 
rors in flux measurements. 

We obtained the electron densities (N c ) from the 
[SHW6717//16731 fine ratio. Although [OII]^3726,29 emission 
lines are brighter than [SII] A67 17,31, the spectral resolution is 
at the limit to resolve accurately the [Oil] doublet lines. [Oil] 
densities are roughly consistent with [SII] densities within the 
errors. 

We derived the electron temperature values of TJOIII] us- 
ing the [OIII]^4363M4959,5007 line ratio. T e [OII] was calcu- 
lated from the relation between [Oil] and [OIII] electron temper- 
atures given by Pilyugin et al. (2006). We were able to measure 
the faint auroral line [OIII]/i4363 for 16 fibers with high S/N 
leading to relative errors lower than 10%. These fibers cover a 
projected area of nearly 19 arcsec 2 equivalent to ~ 0.4 x 0.3 
kpc 2 , including the central starburst region. This area is shown 
in figure [8] superposed on the Ha emission map. T e [OIII] val- 
ues used to derive the ionic abundance ++ were derived from 
emission line fluxes from the high spectral resolution observa- 
tions in order to avoid any possible contamination of /14363 by 
terrestrial Hgl emission. Nevertheless we compared the T f [OIII] 
values derived from both our low and high spectral resolution 
observations, finding that the measurements coincide within the 
errors. 
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Figure [9] shows the distribution of the electron temperature 
TJOIII] (top left panel), electron number density derived from 
the [SII]/16717,31 emission lines (top right panel) and the total 
oxygen abundance spatial distribution (bottom panel). The errors 
are quoted in parenthesis. We show only the fibers where the 
[OIII]/l4363 emission line could be well measured. 

We did not find any statistically significant variance for the 
electron number density; N e values are of the order 100 crrT 3 ap- 
proximately, within the errors. We saw that there are some varia- 
tions in the derived electron temperature (1 1700 - 16000 K). The 
derived values of TJOIII] versus the error in the [ODTJ/14363 
line were plotted and no correlation was found, indicating that 
no significant noise contribution near /14363 was included in the 
flux measurements. 

With regards to the oxygen ionization correction factor 
(ICF), a small fraction of O/H is expected to be in the form 
of the 3+ ion in the high-excitation HII regions when the 
Hell/14686 emission line is detected. We have a measurement of 
the Hell/14686 emission line in 3 fibers. According to the pho- 
toionization models from Stasiriska & Izotov (2003), the 3+ /0 
can be in the order of 1% only in the highest-excitation HII re- 
gions [0 + /(0 + + 2+ ) < 0.1]; therefore, taking our abundance 
results into account, this ICF correction was assumed to be neg- 
ligible in our case. 

The oxygen abundance derived for the 16 positions studied 
does not show correlation with the ionization degree over the 
whole range in 2+ /0 + . There are some variations in the derived 
oxygen abundance in the range 12+log(0/H)= 7.65 to 8.05. The 
error-weighted mean of the derived oxygen abundance values is 
12 + log(0/H) = 7.86 ± 0.05 (lcr). 14 out of 16 O/H measure- 
ments are consistent with this mean value at 2cr (+ 0.1 dex). The 
other 2 measurements present values of O/H consistent, within 
the errors, with + 3<x of the weighted mean. Bearing in mind 
that the errors quoted in the oxygen abundance do not account 
for some observational uncertainties (e.g. pointing errors, seeing 
variations, etc) nor for errors associated to the reddening correc- 
tion and flux calibration, we must assume that the variations in 
the derived oxygen abundance may not be statistically signifi- 
cant. 

In our spectra, neon is seen via the [Nelll] emission line at 
/13869 A. For this ion, we took the [OIII] electron temperature 
as representative of the high excitation zone. We calculated the 
neon abundance assuming that 

Ne/O = Ne 2+ /0 2+ 

Regarding [SII] and [Nil] temperatures we assumed the ap- 
proximation T e [SII] ~ T e [Nil] « T e [OII] as valid, since no auro- 
ral line could be measured in the low excitation zone. The N/O 
abundance was derived under the assumption that 

N/O = N + /0 + . 

Regarding N/O and Ne/O abundance ratios, evaluating the 
contribution of all observational errors to the derivation of 
these abundances, it can be concluded that there is no statisti- 
cal evidence for any systematic variation of N/O and Ne/O [- 
1.15+0.09 < log(N/0) < -0.90±0.11; -1.07±0.18 < log(Ne/0) 
< -0.84+0.10] with O/H for the range of abundances studied in 
this work. The values of the error-weighted mean for log(N/0) 
and for log(Ne/0), -1.04 + 0.06 dex and -0.97 + 0.08 respec- 
tively, are consistent, within the errors, with the log(N/0) and 
log(Ne/0) solar values, -0.86 + 0.12 and -0.82 ± 0.11 (Lodders 
2003), with N/O being slightly (~ 0.1 dex) lower than so- 
lar. The error-weighted mean value for log(N/0) in IIZw70 is 




Fig. 8. [OIII]/15007 flux map. The blue box marks the zone of 
the galaxy where we could measure the [OIII]/l4363 emission 
line. The size of the image is 25" X 21" with a pixel size of 1". 
North is up, east is left. 



significantly higher (~ 0.4 dex) than the mean for metal -poor 
HII galaxies [log(N/0) = -1.46 + 0.14 for galaxies with 7.6 < 
12+log(0/H) < 8.2; Izotov & Thuan 1999]. Still the N/O ratio 
of IIZw70 is consistent with values found for a few HII galaxies 
in the literature (Nava et al. 2006, Izotov & Thuan 1999) with 
similar metallicity. As regards Ne/O abundance ratio, the error- 
weighted mean value for log(Ne/0) derived here (-0.97 + 0.08) 
is lower than the mean value of -0.72 ± 0.07 obtained by Izotov 
& Thuan 1999. However, they also found some galaxies with 
log(Ne/0) similar to the values derived in this paper. In any case 
we should bear in mind that neither the literature long-slit mea- 
surements we are comparing to nor our IFU mosaic cover the 
entire galaxy. 

3.4.1. Physical properties and chemical abundances from 
the integrated spectrum of IIZw70 

In this section some of the more interesting spectroscopic prop- 
erties of the ionized gas in IIZw70, derived by analysing the inte- 
grated spectrum are described. We obtained the integrated spec- 
trum of IIZw70 by summing the emission from each fiber within 
an area of ~ 300 arcsec 2 (~ 1.8 kpc 2 ), enclosing the whole neb- 
ular emission. Integrated spectra were obtained for low and high 
spectral resolution. 

The analysis of the integrated spectrum presents several ad- 
vantages. Contrary to previous studies that attempt to describe 
the average ionization conditions in the ionized nebula using in- 
dividual spectra taken at different locations, we were able to de- 
scribe the average spectroscopic properties over an extended re- 
gion. Another advantage is that the integrated spectrum does not 
suffer from systematic effects that may be found for the spec- 
tra obtained with smaller apertures (Kewley et al. 2005). These 
apertures are usually centered on bright SF knots, presenting 
a bias towards young populations (as a matter of fact in most 
spectroscopic surveys of HII galaxies carried out up to now). 
Furthermore, one can compare the real integrated properties of 
the local SF dwarf galaxies with the properties of SF galaxies in 
intermediate/high redshift for which only their integrated char- 
acteristics are known due to their distance. 

The emission lines corresponding to the integrated spectra 
were measured following the same procedure as in Kehrig et 
al. (2004). Reddening corrected line intensities, normalized to 
H/3 are shown in table [3] as well as the value of C(H/3), esti- 
mated using the Ha/H/3 and Hy/H/3 ratios, the reddening cor- 
rected flux of Ha, F(Ha), EW(H/3), EW(Hff), and corrected ve- 
locity dispersion (see next section for details) derived from H/3 
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Fig. 9. Top left panel: electron temperature distribution (in 
10 4 K) from the [OIII] ^4363/(^4959+/l5007) line ratio. Top 
right panel: electron density distribution (in cm -3 ) from the 
[SII]/16717//16731 line ratio. Bottom panel: the distribution of 
the derived oxygen abundance, 12+log O/H. Note that the re- 
gions shown are outlined by the blue box in figure [8] 

and [OIII]/l4959. The effect of the underlying absorption in the 
hydrogen lines appears not to be important; a 0.3 A EW cor- 
rection for the underlying absorption on the integrated spectrum 
would be enough to bring the observed Balmer ratios to their 
theoretical value within the errors. 

From the reddening corrected line intensities of the inte- 
grated spectrum we derived ionization structure indicators, phys- 
ical conditions and chemical abundances (see table @). These 
quantities were derived as described in section 3.4 for individ- 
ual fibers. 

3.5. Kinematics of the ionized gas 

One of the advantages to using two-dimensional spectroscopy 
is that the spatial distribution of the kinematics of the ionized 
gas can be studied when using high enough spectral resolution. 



Table 3. Reddening corrected line intensities normalized to H/3 
for the integrated spectrum of IIZw70. The reddening coeffi- 
cient, C(Hy3), the reddening corrected flux of Ha and H/?, F(Ha) 
and F(H/3) in units 1CT 13 erg cirT 2 s _1 , equivalent width of Ha 
and H/3, (EW(Ha) and EW(HyS) in A), and the corrected veloc- 
ity dispersions derived from H/3 and [OIII]/14959, (cr cor (H/3) and 
cr eor ([OIII]/l4959) in km s~'), are included. Errors are quoted. 



A Ion flux 



3726 [Oil] 


1.03 ±0.01 


3729 [Oil] 


1.38 ±0.02 


3868 [Nelll] 


0.31 ±0.01 


3889 H8+Hel 


0.19 ±0.01 


3970 He+[NeIII] 


0.08 ±0.01 


4101 H<5 


0.24 ±0.01 


4340 Hy 


0.45 ±0.01 


4363 [OIII] 


0.07 ±0.01 


4471 Hel 


0.06 ±0.01 


4861 H/? 


1.00 ±0.02 


4959 [OIII] 


1.20 ±0.01 


5007 [OIII] 


3.76 ±0.03 


6563 Ho- 


2.80 ±0.03 


6584 [Nil] 


0.17 ±0.01 


6678 Hel 


0.03 ±0.01 


6717 [SII] 


0.29 ±0.01 


6731 [SII] 


0.20 ±0.01 


C(H/3) 


0.30 


F(Ho') 


10.72 ±0.14 


F(H/3) 


3.83 ± 0.05^ 


EW(Hq-)(A) 


-221 ± 18~ 


EW(H/?)(A) 


-44 ± 1 


cr cor {W) (km/s) 


19 ±3 


tr ro ,.([OIII]^4959) (km/s) 


21 + 3 



" Assuming HcifHfl = 2.80 for Case B recombination 
(Storey & Hummer 1995) 



Radial velocity and line widths (cr = FWHM/2.35) were ob- 
tained by fitting a single Gaussian to the emission line. The width 
is corrected for the instrumental profiles and for the thermal mo- 
tions, following the formula: 

2 _ 2 _ 2 _ 2 
®cor ~ ® obs ®ins ®~ther 

where cr cor , cr b s , cr ins and cr t i ier are the corrected width, the 
observed width, the width of the instrumental profile and the 
width due to the thermal motions, respectively. The o- ins is ob- 
tained from the profiles of emission lines of lamps, which is ~ 45 
km s _1 in our high dispersion spectra. To derive the width of the 
thermal profile the electron temperature (r e =14600K) derived 
from the integrated spectrum is used as a representative value, 
finding a value of ~ 3 km s _1 to be applied to the [OIII]/l4959 
emission line. 

In figure QJjjthe maps of radial velocity and of corrected ve- 
locity dispersion, cr cor , resulting from single Gaussian fits to the 
[OIII]/i4959 emission line are shown. 

The corrected velocity dispersion o~ cor map shows the maxi- 
mum values associated to the tails region of the galaxy, with cr cor 
decreasing towards the brightest zones. This map also shows a 
conspicuous minimum close to the position of the minimum ra- 
dial velocity. In addition, the corrected velocity dispersion values 
measured from the integrated spectrum (see table [3} are similar 
to those obtained for the brightest part of the galaxy (see fig- 
ure [TOjl. This is because the velocity dispersion for HII galaxies 
is expected to be dominated by the central (core) component as 
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Table 4. Physical conditions, abundances and ionization struc- 
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shown by Telles et al. (2001). The map of the velocity field ob- 
served presents a radial velocity minimum coincident with the 
locus of the central burst. The difference between this minimum 
and the highest values of radial velocities, found in the tails to 
the southwest, is of ~ 60 km/sQThis blueshifted emission could 
be produced by an outflow from the bright burst of the galaxy. 
We assume the outflow interpretation because of the double-tail 
morphology of the ionized gas and the cr cor increasing towards 
the tails region of the galaxy; whereas for a disk, gas velocity 
dispersion is expected to fall with radius (Binney & Tremaine 
1987, Bertola et al. 1995). 



4. Discussion and conclusions 

4. 1 . Ionized gas and massive stars 

We analyzed integral field spectroscopic observations of the HII 
galaxy IIZw70 in order to investigate the interplay between the 
massive stars and the ionized gas. This galaxy is interacting with 
its nearby companion IIZw71. The contiguous streamer of gas 
between IIZw70 and IIZw71, as well as the SF activity seen 
in both, indicate an ongoing interaction and possible gas inter- 
change between these two galaxies (Cox et al. 2001). We studied 
the spatial distribution of properties of the ionized gas (electron 
temperature and density, chemical abundances, dust extinction, 
excitation, kinematics) from different emission-line maps. These 
maps reveal a central starburst (Cairos et al. 2001, 2001) sur- 
rounded by an extended lower excitation, low surface brightness 
ionized component. From our integrated flux in Ha, corrected 



10 15 20 25 
A RA (orcsec) 



Fig. 10. Kinematics of IIZw70: maps of the corrected velocity 
dispersion, cr([OIII]/14959), and of the radial velocity as mea- 
sured from [OIII]/l4959. The lowest values for both, radial ve- 
locity and velocity dispersion coincide with the central burst lo- 
cus. North is up, east is left. Relative right ascension (RA) and 
declination (DEC) are shown in arcsec. 



by extinction, we derived a star formation rate of SFR(Ha) = 
0.3 M /yr for IIZw70 following Kennicutt (1998); this value is 
in agreement with the value we derived using the flux reported 
in Gil de Paz et al. (2003) from CCD Ha imaging. Kewley 
et al. (2002), using long-slit spectroscopic mapping, derived a 
SFR(Ha) = 0.19 M /yfl Different integration areas and C(Hy6) 
values could be possible reasons for the discrepancy between 
our SFR(Ha) measurement and the value found by Kewley et al. 
(2002). 

For the first time we detected the presence of WR stars in 
IIZw70. For the blue bump, our measurements of the equivalent 
width, EW(WRbump) = -2.2 + 0.6 A, and of the ratio 100 x 
I(WRbump)/I(Hj6) = 4.6 +1.1 (section 3.3) are in good agree- 
ment with the predictions of the evolutionary synthesis models 



Starburst99 (assuming a Salpeter IMF, M llpper = 120 M Q and 
Z=0.004 with high mass loss; Leitherer et al. 1999). The mod- 
els predict EW(WRbump) values around -2 A at ~ 3.2 Myr af- 
ter an instantaneous burst. Corresponding predictions for 100 x 
I(WRbump)/I(Hy6) give values around 3.6 Myr. This fact indi- 
cates that the central burst in IIZw70 is very young. Assuming 
an age for the burst of 3.2 Myr (as deduced from our measure- 
ments), we obtained a ratio of WR/O = 0.056 from Starburst99. 
The presence of a bump at ~ 5800 A could be an indicator of 
the existence of WC stars. A broad feature at this wavelength 



3 Cox et al. (2001) found a value of 67 km/s for the rotational speed 
of the HI neutral gas. 



4 We calculated this value converting the SFR(Ha) given in Kewley 
et al. (2002) from the distance they adopted to the distance adopted in 
this paper (see table 1). 
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has been marginally detected in the same spectra, and though 
the measured EW is consistent with the same model predictions, 
its statistical significance is poor due to the low S/N. Overall, a 
total number of ~ 17 WN+WC stars are predicted by the models 
from our findings. 

The nebular emission of Hell/14686 detected in this work is 
shown to be spatially extended, though it is not coincident with 
the location of the WR bumps by a few arcsecjfl This nebular 
Hell emission is seen to the southwest of the box encompassing 
the probable ionizing clusters (see figure 4), suggesting the ex- 
istence of very hard ionizing radiation which, in principle, could 
not be necessarily related to the WR stars located towards the 
northeast corner. 

In the overall picture, IIZw70 is dominated by the presence 
of a young central starburst with a significant SFR. The peak of 
the Ha emission is located close to the probable youngest stellar 
clusters (see figure 0J. These young starburst could be produc- 
ing strong winds which have been able to develop an ionized gas 
outflow from which we see the blueshifted component as a fea- 
ture in the radial velocity map of the galaxy (see figure 10). This 
outflow is thought to have been responsible for the minimum in 
optical extinction (slightly displaced from Ha emission peak by 
~2" to the west) apparent in the Ha/Hp map, where dust could 
have been swept up by the wind. The characteristic size (~ 150 
pc) and velocity (~ 60 km/s) of the outflow zone are consistent 
with previous findings in the literature (see e.g. Roy et al. 1991) 
for the giant HII complex NGC 2363). We can also see that the 
area with lower dispersion velocity partially overlaps with the 
zone of the gas outflow, between the probable ionizing clusters 
location and the zone of minimum extinction. From the kinemat- 
ical point of view, this picture appears to be consistent with the 
velocity dispersion map. The narrowest line profiles are domi- 
nated by the bulk of the emission from the central HII regions, 
as expected (Telles et al. 2001), whereas the rest of the ionized 
gas, with a lower surface brightness, seems to be suffering the ef- 
fects of the winds and its kinematics should be associated to the 
structures generated by the likely interaction of the winds with 
the ambient medium. This situation is reminiscent of what has 
been found for other extragalactic HII regions and HII galaxies 
(e.g. Munoz-Tunon et al. 1996, Telles et al. 2001) which host 
strong starbursts. 

4.2. Chemical enrichment and abundance variations 

In this section, based on our chemical abundance measurements, 
we present an analysis on the chemical enrichment of the ionized 
gas and chemical abundance variations in IIZw70. 

The oxygen abundance derived here for the integrated spec- 
trum of IIZw70 is 12 + log(0/H) = 7.83 + 0.04. Previous works 
reported direct values of oxygen abundance for IIZw70 using 
long-slit spectra positioned in different places across the bright- 
est regions of IIZw70. Lequeux et al. (1979) quoted a value of 
8.07; Kobulnicky & Skillmann (1996) derived an oxygen abun- 
dance of 8.07 + 0.08 whereas Shi et al. (2005) found 12+log 
(O/H) = 7.69. In principle, the differences between our O/H 
measurement and the values found in the literature could be ex- 
plained by different integration areas (see figure 9, bottom panel; 
Kewley et al. 2005) and atomic parameters used. 

Thus far, in most works on chemical abundances in dwarf 
galaxies there is no clear evidence of abundance variations. 
Though selected examples of nitrogen enhancement in dwarf 

5 We should bear in mind, however, that a small contribution of neb- 
ular HeII/i4686 could have remained unnoticed in the broad feature 



galaxies have been reported in the literature, notably the case 
of NGC 5253 (e.g. Lopez-Sanchez et al. 2006 and references 
therein), this is not the case of O/H. Many works in the literature 
have searched for a possible oxygen abundance variation within 
dwarf galaxies. Kobulnicky & Skillman (1996) presented optical 
spectroscopic chemical measurements of the ISM in the irregular 
galaxy NGC 1569 that reveal no substantial localized chemical 
self-enrichment (i.e., "local pollution"; see e.g. Kunth & Sargent 
1986; Pagel, Terlevich & Melnick 1986; Pilyugin 1992). Vflchez 
& Iglesias-Paramo (1998), using bi-dimensional, long-slit spec- 
troscopy, showed that IZwl8 presents a substantially homoge- 
neous chemical composition over the whole galaxy. From long- 
slit spectra of 16 HII regions in the dwarf irregular NGC 1705, 
Lee & Skillman (2004) reported that there is no significant spa- 
tial variation of oxygen abundances for this galaxy. Izotov et al. 
(2006), from two-dimensional spectroscopy of the HII galaxy 
SBS0335-052E, found a variation of 0.4 dex in the oxygen abun- 
dance. They concluded that these variations may not be statisti- 
cally significant due to unaccounted error estimates. Lee et al. 
(2006) obtained a slightly significant (3.2cr) oxygen abundance 
gradient of -0. 16 + 0.05 dex kpc 1 within the Local Group dwarf 
irregular galaxy NGC 6822. However they claim further deep 
high-quality spectra of nebulae and stars are needed to distin- 
guish clearly between either a zero or a nonzero slope. 

In this work we calculated the O/H abundance in 16 positions 
across IIZw70 directly from the measurement of [OIII]/14363. 
From these values we found a maximum variation in the derived 
oxygen abundances of ~ 0.4 dex, similar to the range reported 
by Izotov et al. (2006) for SBS0335-052E. In order to study 
whether this variance is statistically significant for an effective 
abundance variation we must include all other unaccounted un- 
certainties (variable seeing, flux calibration, etc), which surely 
contribute to the total oxygen abundance error but are difficult to 
estimate. In any case, our observations indicate that, within ~ + 
0.2 dex, the ionized gas in the central burst of IIZw70 appears to 
be chemically homogeneous in O/H over spatial scales of hun- 
dreds of parsecs. In the case of the neon-to-oxygen ratio, Ne/O, 
the range of values as derived from our point by point abundance 
analysis, does not show inhomogeneity within the errors; as is 
also the case of the nitrogen-to-oxygen ratio, N/O (see section 
3.4). All these findings give us the upper limits to any chemical 
inhomogeneity in the ionized gas in the brightest part of IIZw70. 
In accordance with these results we can say that we are probably 
not detecting a "local pollution" case in IIZw70, consistent with 
most results obtained by other works cited previously. 

A favourable explanation to this apparent degree of homo- 
geneity in dwarf galaxies is given by the scenario presented in 
Tenorio-Tagle (1996). According to this scenario, in those places 
in which one expects that newly produced metals could be re- 
turned to the ISM, the "local pollution" is not detected in the 
warm phase of the ISM. Rather the newly synthesized metals 
are returned into the hot phase of the ISM, and there is a sig- 
nificant time delay before these newly produced metals can be 
detected in the warm phase of the ISM. Variations of the oxy- 
gen abundance in the ISM of HII galaxies could give constraints 
to the physical mechanisms involved in the recycling processes 
(e.g. Recchi et al. 2001). 

Analysing the predictions of chemical evolution models, we 
verified that continuous SF models (Molla & Diaz 2005) pre- 
dict an enhancement in O/H of up to 0.4 dex, to be produced 
in spatial scales of ~ 1 kpc along their disk, over timescales 
of 8 Gyr. The important point here is that these models predict 
O/H enhancements that could be achieved keeping the N/O ra- 
tio approximately constant throughout the process, and at an ab- 
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solute value consistent with our observations (see section 3.4) 
to within 0.1 dex. Further observations with IFUs in 10 m class 
telescopes and additional measurements of chemical abundances 
in IIZw70, covering a larger area, should be required in order to 
draw a firmer conclusion on the statistical contribution of local 
variations to the chemical homogeneity. Given the interest in the 
problem of chemical homogeneity of dwarfs galaxies and its role 
in their chemical evolution (e.g. Legrand et al. 2001, Recchi et 
al. 2001) a clearer view of the observational situation is probably 
warranted. 
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